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Using a theory of polarizable fluids, we extend a variational treatment of an excess 
electron to the many-electron case corresponding to finite metal concentrations in 
metal-ammonia solutions (MAS). We evaluate dielectric, optical, and thermodynam- 
ical properties of MAS at low metal concentrations. Our semi-analytical calculations 
based on a mean-spherical approximation correlate well with the experimental data 
on the concentration and temperature dependencies of the dielectric constant and 
the optical absorption spectrum. The properties are found to be mainly determined 
by the induced dipolar interactions between localized solvated electrons, which result 
in the two main effects: the dispersion attractions between the electrons and a sharp 
increase in the static dielectric constant of the solution. The first effect creates a 
classical phase separation for the light alkali metal solutes (Li, Na, K) below a criti- 
cal temperature. The second effect leads to a dielectric instability, i.e., polarization 
catastrophe, which is the onset of metallization. The locus of the calculated critical 
concentrations is in a good agreement with the experimental phase diagram of Na- 
NH3 solutions. The proposed mechanism of the metal-nonmetal transition is quite 
general and may occur in systems involving self-trapped quantum quasiparticles. 
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I. INTRODUCTION 



Metal-ammonia solutions (MAS) have been the subject of numerous experimental and 
theoretical studies, since their discovery by Weyl in 1864. Many experiments have been 
performed and a large volume of experimental data have been accumulated (for a review, 
see [ll, y, y, Q] ) . However MAS still remain essentially an unsolved problem. The interest in 
this fascinating system is twofold: first, it exhibits a rich variety of different phases includ- 
ing a phase separation and a metal-to-insulator transition, and second, the the electronic 



concentration is low in the metallic state so that correlation effects are important. 
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'his 



intriguing system therefore remains of sustained interest sl, 0, ^, ol, Q, 

Once an alkali metal is dissolved in liquid ammonia, it immediately dissociates to give 
separated entities with unlike charges: the solvated ions and electrons. At these concentra- 
tions, the solvated electrons form a cavities free of solvent due to short-range interactions 
with ammonia molecules, in which they localize with the help of the polarization carried 
by the surrounding ammonia molecules [l^. Path integral simulations Q, 0, Q, Q and 
density functional approaches 2^, 2l| , providing an evaluation of the microscopic structure 
around solvated electrons, indicate also the possibility of solvated dielectrons, i.e. paired 
electrons trapped in a cavity, eventually as metastable states. The solution remains non- 
metallic (electrolytic) at these low concentrations and has an intense blue color independent 
of the type of alkali metal, while the optical absorption spectrum does not change up to 
concentrations of about 0.1 mole percent of metal (MPM) 221]. For metal concentrations 
exceeding 7 MPM, the MAS is a liquid metal with a typical bronze coloration. The con- 
ductivity of MAS reaches that of liquid mercury at concentrations of about 20 MPM [l|. 
However, the blue and the bronze phases are macroscopically separated in the intermediate 
range of concentrations varying from 1 to 10 MPM, resulting in a miscibility gap below a 



critical temperature 



23 



24j . Numerous experiments give evidences about a metal-nonmetal 



(MNM) transition in this range (also called metal- insulator transition (MIT) in the paper), 
indicated by the evolution of the conductivity, the Hall resistance, the Hall mobility, and 
the Knight =Mft fl Q. h. part.cula, the conductivity of Na-NH, .o.ut.on. inc.ea=e= by 
three orders of magnitude between 3 and 6 MPM pj. 



Previous models have considered the Mott mechanism 



25j or have involved an association 



of localized electrons in clusters 26|] to explain the MNM transition, but were not able to 
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explain satisfactorily the whole phase diagram. Concerning the Mott scenario, for example, 
which is the usual MIT mechanism in doped semiconductors, the MAS seem to be an 
exception because a phase separation arises in the concentration range where the MNM 
transition occurs. We have recently shown that the MNM transition is rather driven 
by the old Goldhammer-Herzfeld (GH) mechanism of metallization 28|, l29|]: the quantum 
fluctuations of the electrons localized in cavities induce a polarization catastrophe resulting 
in metallization of the solution. This proposal follows our previous studies of polarons 



interacting at low densities, for which a similar scenario may occur 30|, l3l|, l32| , and even an 
analogy between high- Tc cuprates and MAS was also proposed 

We emphasized the prominent role of the dipolar-dipolar interactions i.e. van der 
Waals interactions in the sense of the London dispersion forces in the MNM transition. 
Initially, the idea seems unlikely, since these forces are usually very weak with respect to the 
Coulomb interactions in ionic systems. In fact, the situation is more complicated. Diluted 
MAS are simply electrolytes consisting of solvated cations (metal ions) and solvated electrons 
playing the role of anions in the solution. The two species have a different behavior, (classical 
or quantum), depending on which energy (or time) scale is considered. Both species behave 
classically at long time scale. The same is true of the degrees of freedom: the cation 
coordinates and the center of mass coordinates of each solvated electron may be treated 

n 

classically. We note [27] that at the experimental densities of the NMM transition (~ 4MPM) 
the Debye length is estimated to be lA that is lower than the mean distance between the 
ions (~ I2A) by an order of magnitude. The electrostatic Coulomb interactions are thus 
considerably screened and may not be responsible of the MNM transition. Overlaps between 
the localized electron wave functions also remain negligible at this concentration. Then the 
next step is to consider polarizabilities of the different species: solvent, cations, and solvated 
electrons, which are respectively given by: aNH^ ~ 18.8 aj] [sj], ai^a+ ~ 1-34 [sS^ (for 
sodium metal), ao{0) = / mujQ{TY ~ 913 aj] (in Born units ao = h'^/me^, where m and 
e are respectively the electron mass and charge). The huge polarizability of the solvated 
electrons is due to the internal electronic transition in the cavity. It occurs roughly around 
Uq ~ 0.9 eV, leading to the above polarizability. Notice that these electronic transitions are 
of Franck-Condon type since the relevant time scale Uq^ is much shorter than the relaxation 
time T of the solvent around localized electrons, which is about several picoseconds 36|. 



This polarizability of the solvated electrons induces instantaneous forces of two kinds. First, 
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there are dip ole- charge forces: the interaction between the induced dipole momentum of 
electrons and the classical degrees of freedom. These forces vanish on average due to the 
spherically symmetric distribution of charged species. Secondly, it also induces interaction 
between quantum fluctuating dipoles, i.e. dispersion interactions. Since the pioneering work 
of Onsager, it is well-known that these forces are not additive and must be collectively 
treated. This theoretical task have been achieved by several groups and formulated as the 



theory of quantum polarizable fluids 



37 



38, 



39, 



40 



41 



42 



43|, 



44| . We extensively use these 



results in order to evaluate the contribution of the dispersion forces between the solvated 
electrons and their physical consequences. Mainly, the dipole-dipole interactions are self- 
induced cooperatively, in contrast to to the case of Coulomb forces which always tend to 
be self-screened. That is why the dipolar forces may induce a polarization catastrophe at a 
finite concentration of metal, provoking the onset of metallization. This scenario essentially 
follows the GH approach except that it occurs in a system having charges: the solvated ions 
and electrons. 

Applying a semi-phenomenological treatment, we propose the excess free energy to con- 
tain three kinds of terms, namely, an additive contribution of noninteracting electrons, the 
excess free energy due to the classical degrees of freedom, and the van der Waals part due 
to the dispersion interactions. In the paper, we give a complete theory of the nonmetal- 
lic phase, including details which were beyond the scope of our previous publication {2^. 
We will evaluate thermodynamic quantities, dielectric, and spectroscopic anomalies of the 
solution. Although the complete theory of the phase diagram requires the consideration 
of coexistence of the metallic and the nonmetallic phases, we restrict ourselves here to the 
nonmetallic phase and focus on the properties of excess electrons in this phase. The metal- 
lic states and a crossover regime will be considered in our next papers 45|]. The layout 
of this paper is the following. In Section [TTl we describe our model based on a variational 
semi-continuum treatment of MAS and accounting nonpolar, electrostatic, and dipolar inter- 
actions between localized electrons. In Section UTTl we present the results obtained within the 
framework of our model. The last section is devoted to discussions, future extensions, and 
applications. Two appendices contains special topics: derivation of the effective free energy 
functional and evaluations of the effective polarizability for interacting electrons. Atomic 
units are used throughout. 
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II. MODEL AND APPROXIMATIONS 



A. General outline of the problem 



Upon dissolving the metal, a broad optical absorption appears around loq ~ 0.9 eV, 
which does not depend on nature of the solute. This experimental fact indicates that solvated 
electrons dissociated from metal ions are responsible for the absorption and form independent 
entities |l| . The problem of a single solvated electron has been intensively studied (see Sec. 
Ill Bl) . however the many-body problem posed by a finite concentration of solvated electrons 
is a more difficult task, because it includes classical and quantum correlations, as well as 
short and long-range interactions. Indeed, there are few theoretical approaches to investigate 
electrons solvated in MAS at finite metal concentrations except the ones excee ding 10 MPM, 
where the system is highly metallic and localized electrons have disappeared 46|. 

We propose a semi-phenomenological model following which we treat the main relevant 
physical features of the diluted MAS (sketched on Fig. [1]). As will be discussed later (Sec. 
Ill Bl) . the single solvated electron is localized in a cavity of diameter a and surrounded by 
solvent molecules. The number of molecules which compose the solvated cavity is estimated 
to be about of 8 in diluted solutions (< 5 MPM) l|. The individual cavity diffuses as a whole 
in the solution due to the thermal fiuctuations exactly as would do a classical negatively 
charged ion. The variables {rI^-'^} determine their cartesian coordinates, while {R^^-*^} are 
the cartesian coordinates of the solvated cations. For sake of simplicity, we assume that 
each solvated cation and each solvated electron have the same diameter a. Thus, from an 
electrostatic point of view, the diluted MAS behaves as an electrolyte consisting of cations 
and anions with diameters a dissolved in a solvent with dielectric constant e{uj = 0) = eg, 
which is the one of pure liquid ammonia. Such a solution may be statistically treated by 
the restricted primitive model (RPM), while an analytical expression of the free energy can 
be found by within the RPM (see Sec. HTCl) . 

But this is not a complete description (see Appendix 1). In addition to this classical part, 
there are extra degrees of freedom for the solvated electrons due to the possible individual 
electronic transitions inside their solvated cavities. We introduce variables {u^^J'} which 
are the relative cartesian coordinates of the electrons with respect to their centers of mass 
{Ri^^}. The de grees {u''^^} are quantum mechanical in contrast to the variables {Ri.^-*^}. 
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This separation between the quantum and the classical degrees of freedom may be viewed 
as a Born-Oppenheimer approximation justified by the relevant energy scales, i.e. fc^T for 
the classical degrees and Huq for the quantum ones, which magnitudes roughly differ by 
two orders. Physically, it means that the duration of an electron transition from its local 
ground-state to an excited state (~ uJq^) is sufficiently small with respect to the relaxation 
time of the solvent. Consequently, the quantum correlations or interactions induced by these 
electronic transitions are screened by the high-frequency dielectric constant of the solvent, 
i.e. e{uj ~ oo) = Cqo- This difference of screening for the different interactions, ~ 20 for 
the classical ones, and too ~ 1-7 for the quantum ones, plays a significant role in our the ory 



as it was already identified for the case of the melting of a polaron Wigner crystal 30 



31 



321]. We carry out in Appendix 1 an expansion of Coulomb interactions and show that it 
includes the dipole-dipole interactions between the solvated electrons apart from the classical 
electrostatic interactions. This is nothing but their dispersion interactions as already said 
in the Introduction. They play a dominant role in the thermodynamical behavior of diluted 
MAS (see Sec. IIII CI) owing to the conjugated effects of the large polarizability of solvated 
electrons and the weak screening by the solvent (too as we have seen). Therefore at metal 
density n < 5MPM, the change Afnm{i^) in the total excess free energy per solvated electron 
may be written as the sum of three terms: 

Afnmin) = fo + AUn) + A/,H, (1) 

where /o is excess free energy of a single electron, Afci{n) is the contribution of the classical 
degrees of freedom, and Afa{n) is the van der Waals contribution due to dispersion forces 
between solvated electrons. 



B. Free energy of a single solvated electron 

Numerous theoretical and experimental studies have indicated that 

the polarization due to the surrounding molecules and the cavity formation due to short 
range repulsion plays the main role in the electron solvation. Following this concept, the 
free energy of an excess electron in an infinitely diluted MAS may be written as 
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where the first term Pg/2 is the quantum mechanical kinetic energy of localization for the 
solvated electron, is the mean momentum of the electron which depends on the cavity 
radius rg. The Heisenberg inequality leads to pe ~ Ck/re, where Ck is a constant of propor- 
tionality which depends on details of the electron density distribution. The second term in 
([2]) is the potential contribution due to the polarization of solvent. This polarization has two 
contributions due to the existence of two sources of polarization (as in the polaron problem). 
The first depends on the electronic polarizability of the solvent molecules (ammonia) and 
provides the high frequency dielectric response Eoo- The second corresponds to the orienta- 
tional polarization of the ammonia molecules. It provides, together with the first, the low 
frequency (static) dielectric response Es of the solvent. The last two phenomenological terms 
in ([2]) are the nonpolar contributions due to the cavity formation, i.e. the work necessary to 
perform the cavity formation. It contains a surface term (A is the surface tension) and a vol- 
ume contribution corresponding to the work of an external pressure Ps oc (3~^ns. It depends 
on the solvent density n^, temperature and details of the local structure around the 
cavity formation (the size of the first coordination shell, the number and the orientation of 
solvent molecules in the shell, etc). We approximate the influence of all the details by the two 



phenomenon 
theories 48 



ogical parameters Cs and Cy, which can be evaluated within the microscopic 



49( 1 . or adjust by fitting experimental data. The extremum of the free-energy 



/o(re) provides consistent evaluations of the cavity radius r^, i.e. dfn{rp)/drp = 0. Equation 
([2]) is a modified variational treatment of the one-polaron problem 50|] within the framework 
of the semi- continuum approach [isj ]. 

To take into account the internal electronic transition occurring at frequency ujo{T), we 
relate phenomenologically this frequency to the cavity size as 

-o(T) = (3) 

e 

where Cr(T) is a parameter depending both on the temperature and on the local potential 
felt by the internal degrees {uj}. It is well-known that C.,. = 3/2 for a three-dimensional 
parabolic potential, however the constant deviates from this value for anharmonic oscillators. 
The simple analysis of experimental data yields Cr ~ 1.25 for the infinitely diluted MAS at 
T = —70^C [51], taking = 3.2 A and ujq = 0.9 eV. This indicates weak deviations of the 
actual potential from the harmonic oscillator. Hence, the constant C,r(T) can be fitted to 
provide relations between the calculated values of rg and the experimental data on ujo{T). 
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C. Classical contribution to the excess free energy 

The classical part of the excess free energy is composed by charged hard spheres of same 
diameter a dissolved in a solvent with dielectric constant eg. It yields two contributions: the 
nonpolar (A/„(n)) and the coulomb (A/c(n)) parts: 

A/eKn) = A/n(n) + A/e(n). (4) 

Analytical expressions for these contributions are well-known. Following the Carnahan- 



Starling result 



5d, 



53|, the nonpolar part is expressed as: 

rj{A — 3r]) 



ln(n/?2s) — 1 + 



(5) 



where is the solvent density, while rj = irna^/S is the packing factor. In the above 
expression, we take into account that the density of hard spheres is equal to 2n (cations 
plus solvated electrons). The first term is the ideal part, while the last term is the hard- 
sphere exclusion part. Since for the densities that we consider r/ << 1, the ideal term is 
preponderant. For example, at n = 4 MPM, t] ~ 0.04, and only few terms in the low-density 
expansion of are physically relevant. 

For the electrostatic contribution, we use the results of the mean spherical approximation 
(MSA). The contribution Afc{n) is expressed in terms of the inverse screening length 7 = 
([1 + 2(t/£d]1/2 _ iy2a, where the Debye length is in = {Snpne^ /es)'^^^. The Coulomb free 



energy part is then given by 54l |: 



/?A/cH = - r 1^ -11 + v^- (6) 

D. Van der Waals contribution to the excess free energy 

The most difficult part of the free energy to evaluate, comes from the contribution of the 
dispersion forces between solvated electrons. Because we have separated the quantum and 
the classical variables, we may write the quantum part as the hamiltonian of interacting 
Drude oscillators (Appendix 1). Introducing the dipolar momentum nij = euj, the quantum 
part of the Hamiltonian for a given configuration {rI^-^} is written as: 

N 2„2 1 1 N 

H = Y. + ^ + ^E-^^(|R- - Ri|) + ^Em. ■ T(RL - Ri) ■ m,, (7) 
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where tTj is the conjugate of the variable rrij, and = /mujQ is the bare polarizabihty 
of a solvated electron, Uhsir) is the hard-sphere interaction potential Uhsij" > cr) = and 
Uhsir < o") = oo, while T(r) = 3r~^rr — r~^\ is the dipolar rank-2 tensor. 

To evaluate different physical properties, two quantities are important: the density of 
states per particle D(ct;, {r1^^}), and the average polarizabihty a{uj, {rI^^^}). Both are 
related to the eigenvalues {co'A({Ri^"'^})} of the collective modes of ([7]) and depend on the 
given configuration {rI^-*^}: 

d(..,{r1^^}) = ^E^(^-^a({r1^^})), 

Using mathematical relations of the theory of distributions, namely, /m(l/(u; — zO)) = tt5{uj) 
and dioj"^ — oJq) = 6{uj — ujq)/2uj, we can directly relate the above quantities: 

D(.;,{Rl^>}) = '^"-^"'^'^-'''» , (9) 

71 

where ai{uj, {R^^}) is the imaginary part of a{uj, {r1^^}). 

The next step is to proceed the thermodynamical average over the classical degrees of 
freedom in order to define the DOS and the effective polarizabihty of a solvated electron: 

Dicj) = (d{u;, {rI^^^})) , aH = (a{u;, {rI^^})) , (10) 



where the brackets (■ ■ ■) denotes the thermodynamical average over the variables {rI^^}. 
Rigourously, the thermodynamical average should also include average over the positive ions 
({R+^^}), taking into account the correlations between classical particles due to the Coulomb 
interactions. However, as we have discussed in the Introduction, the Debye length is very 
small (~ lA ) at the densities that we are considering, so that we neglect the correlations 
due to the Coulomb interactions in the thermodynamical average, taking only into account 
the correlations due to the dipolar and the nonpolar interactions. Under this assumption, 
we are left with a polarizable neutral fluid represented by ([7]) , in which the classical variables 
now implicitly correspond to the positions of neutral classical hard spheres. 
In the framework of the theory of polarizable fluids, the problem has been extensively 



studied 



37 



38 



39 



40 



41 



42 



43l . l44j ]. we use this approach and refer the reader to these 
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references for more details. The main result is that the effective polarizability a{uj) defined 
by ffTOj) is calculated via the self-consistent equation: 

(oo) = «o"' (^) - 2f/(« (io) /eoo)/eoo, (11) 

where the last term accounts for the correlations between induced dipoles,while aoi^uj) ~ 
(ti^o ~" 1^^)"^ is the bare frequency polarizability of the particles. Equation ( fTTl) has been 



derived in [39|] , we only modify it by taking into account the high-frequency screening by the 
solvent (the use of eoo in (iTTl) ). The function Sail (a) / (3 corresponds to the dipolar part of the 
internal energy per particle of a classical liquid of nonpolarizable particles with permanent 
dipole momentum (3a;//3)^/^ 37|,l39|. Once this function is known, all the physical properties 



of the system can be evaluated by solution (fTTj) . The simplest method to obtain U{a.) is the 



Pade approximation 39|, |40||, which is an interpolation between the case of low and large 



polarizabilities a{ijj) (Appendix 2). Adapting this method to our case, we get: 

TT( ( \l \ h{nrl)na{uj) 

where Io{x) and Ii{x) are analytical functions depending on dimensionless density x = Snr^ 
[s^]. As a result, we obtain the quadratic algebraic equation for a{uj): 

[eoo/i(a;) - 2Io{x)nao (u;)]a^ (uj) = [8e^rf - eooIi{x)ao{Lj)]a (u) - 8e^rfao(tu), (13) 

whose the complex solution is: 

a{uj) = ar{uj) + iai{u). (14) 

The imaginary part ai{uj) is nonzero only in a finite range of frequency uj^(T,n) < uj < 
u;+(T, n). This is a direct consequence of the Born-Oppenheimer approximation made in 
the model (separation between classical and quantum degrees of freedom, as it is also done 
in the theory of polarizable fluids). In general, all quantities may be evaluated numerically 
by the above equations. However, all the quantities may be found analytically in the special 
case of dilute solutions for which nrg <^ 1 (see Appendix 2). For example, the imaginary 
part of the effective polarizability is given by: 

= — , (15) 
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where Lj±{n) are the edge frequencies: 



In that case the DOS {D{uj) ~ ujai{uj)) has a semielhptic form. 

Once the above calculations have been performed, it is easy to evaluate the van der Waals 
contribution to ([1]), which is directly given by the difference between the free energy of the 
interacting Drude oscillators and the free energy of the same but non-interacting oscillators 



55l | . For the temperature of interest T < 250 K and the bare frequency of the dipoles c^o ~ 1 



eV, we have jSuo ^ 1. Then, the free energy per oscillator in the non- interacting case is 
just the zero-point energy 3huJo/2. The free energy is related to the DOS so that the van 
der Waals contribution is just given by: 

A/. = i ln(2 snM^])d^ - (17) 

This equation must be evaluated numerically together with a{uj) and D{uj). However in the 
diluted solutions, we may obtain an analytic expression (Appendix 2): 

Af, = -u;o[^^,+0{n^atr:% (18) 

This contribution is negative as it should be and the total free energy Fd ~ nAfd is propor- 
tional to n^, as in the case of the van der Waals model for the liquid-gas transition. 



E. Dielectric and thermodynamical properties 

To make the link with the dielectric properties, we calculate the total dielectric con- 
stant of the solution as the one of a set of independent spherical particles with effec- 
tive polarizability a{uj) in a medium with dielectric constant eNHa{T,uj) (pure ammo- 
nia). The total polarization of the sample under the macroscopic field E{uj) is by def- 
inition Piot(cu) = [(e(T, tu) — l)/47r]E(co') where e(T,u) is the requested dielectric con- 
stant of the solution. The polarization may be written as the sum of two contributions: 
Ptot(t^) = PNHai^^) + AP{u) due to the solvent and the polarizable spheres, respectively. 
If we assume that the dielectric properties of the solvent in itself are not affected by the 
presence of the polarizable spheres, we may write P^Hsi^) = li^NHsiTju) — l)/47r]E(c<j), 
while the dipolar contribution to the polarization is just given by: 

AP(^) = (£(^li^^l^WT^E(..). (19) 
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At the same time, the excess polarization due to the quantum polarizable spheres in the 
solvent is expressed as (see Appendix 2): 

AP(cu) = n«(u;)E,(cu)/eoo, (20) 

where we have assumed that only the high dielectric constant of the solvent screens the 
polarizability of the spheres (due to their quantum nature), and E£{uj) is the local field 
viewed by the given set of polarizable spheres. The Lorentz local field expression gives: 

Ei{uj) = E{uj) H , Ei{uj) = E{u) H . (21) 



Using (HM . (I2U]) . and ([21]), we obtain the modified Maxwell-Garnett expression |56|): 

e{T,uj) - eNHaiT,^;) 47r 



e(T, uj) + 2eNHi i^, uj) Se^o 
which gives by inversion the expression of the dielectric constant of the solution as function 
of its constituents, i.e. the solvent and the solvated electrons. Notice that in the expression 
of the total polarization we have ignored the polarization due to the presence of the classical 
mobile charges, since it is not vanishing only at very low frequencies due to the very low 
mobility of the ions in the electrolyte. Our expression (!22l) is thus valid provided to > Tj^^ 



where typically r^) ~ 10 s 



Concerning the thermodynamic properties, we may calculate the excess pressure, the 
excess chemical potential, and the excess compressibility k due to the presence of the solvated 
electrons in the solutions by the usual relations: 

Ap = n —— , A/i = Afnm + n — , k = n——. (23) 

on on on 

To complete the calculations, we also mention that the solvated radius Ve may be self- 
consistently obtained by minimizing the excess free energy ([1]): 

= 0. (24) 

F. The use of experimental data to fit plienomenological parameters. 

Our input phenomenological parameters are C5, Cy, C^, CriT), a, eoo, and egiT). The 
first four constants can be derived from the experimental data on the absorption spectrum 
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z/(a;) and on the excess molar volume directly related to r^. The mean momentum is also 
related with the first moment of the absorption spectrum, i.e., Pe = uji = J ujv{uj)duj. 

In general, the dielectric and thermodynamical properties of the localized electrons should 
be calculated self-consistently, but it leads to complicated concentration and temperature 
dependencies of the cavity size re{n) and of the maximum uJo{n) of the absorption spectrum. 
To avoid these cumbersome calculations, we fix these quantities re(n) = r^. (obtained for 
n ^ 0) and ujQ{n) = ujq (also obtained for n ^ 0) at the first step of our calculations. 



Using experimental data on //(cu) and Tg at zero metal concentration and T = — 70*^C 51| . 
we obtain Ck ~ 1.5 and Cr ~ 1.25, while we set = 1 and Cy = 1-75 to obtain the 
experimental value rg = 3.2 A and uJo{T = — 70°C) = 0.9 eV. 

We also use the experimental data on co'o(T) 58|]. This last parameter leads to an implicit 
temperature dependence of a(uj), we take duoiT) / dT = —2.2-10"'^ eV/K in the temperature 
range {—70°C < T < +70°C) [58| and use the appropriate data on dielectric constants of 
hquid ammonia, i.e., eoo = 1-756, e,(r = -70°C) = 25, and des/dT = -OJ^K"^ [591 . 
Finally we also need data on the density of pure ammonia, deriving them from 
unhs = 0.0255A-3 at T = -70°C, while dnNuJdT = -5 • 10"^ A-^K-\ 



60|, namely, 



III. RESULTS AND ANALYSIS 
A. Criticality of the nonmetallic phase 

Several transitions can arise in nonmetallic MAS due to a variety of classical and quantum 
effects. First, we consider quantum instabilities: one of them is related to the GH scenario 
of metallization 29|], following which Herzfeld took the point of divergency of the dielectric 
constant as the onset of metallization. Such an instability also occurs for the same reason in 
our model, and we call it as n^{T). It depends on T due to the temperature dependence of 
the bare polarizability through the parameter ujo{T). Another quantum instability arises due 
to softening of the collectives modes, that induce a quantum transition at a critical density 
riciiT) also depending on temperature. We will see in Sec. IIIIBI that both instabilities 
riooiT) and nci(T) physically coincide and correspond to the onset of metallization. 

Another important consequence of the model is a thermodynamical instability. It is 
well-known for a classical fluid that the van der Waals interactions compete with the short- 
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range interactions to drive a liquid-gas phase separation below a critical temperature T^^ . 
The same transition occurs in our model. Looking at the locus of this "liquid-gas" phase- 
separation in its low-density regime, the spinodal curve UgiT) is evaluated and discussed 
in Sec. IIII Dl by examining ([T]). The term "liquid-gas" refers, in our model, to a phase 
separation between a phase with a low density of solvated electrons, and a phase with a 
higher density of the same objects. In fact, this classical phase separation is hidden by the 
quantum MIT at rid > n^siT) and should be replaced by a phase separation between the 
insulating phase of solvated electrons and the metallic phase of free electrons. This point has 
)een already partly discussed by us [27|, details will be reported in our future publications 



45|. We only consider here the low density part of the spinodal curve ns{T). 



B. Density of states and stability of the insulating phase 

The effective polarizability a{uj) and the DOS D{uj) can be immediately deduced by ( l37|l 
and ( JTOl) for the known (and fixed) cavity size rg = 3.2 A . The typical behavior of D{uj) is 
shown in Fig. [2]at various concentrations and T = — 70°C, while the dimensionless imaginary 
part ai^uj) / ao{0) and real part ar{uj) / ao{0) of the effective polarizability are shown in Fig. [3] 
under the same conditions. As it is seen, for a given temperature, the DOS and the effective 
polarizability broaden progressively as the concentration n rises. Notice that the maximum 
of the imaginary part of a{uj) shifts to low frequencies as n increases, while the maximum of 
the DOS remains located at the same position. At low enough concentrations the deviations 
of the DOS and the imaginary part of a{u!) from a semielliptic forms are minor, but the edge 
frequencies uj±, given by (fT6|) . monotonically vary with n (see Fig. H]). The DOS broadening 
is roughly proportional to Auj± = uj+ — lu- ^ rn}!'^ at low concentrations. More precisely, 
the low edge lO-in) progressively decreases as n increases, indicating a softening of some 
collective modes of the solvated electrons. This softening with increasing density is directly 



related to the dipolar interactions 30|, l3l|, |6l| . Then, the natural condition of stability is to 
be ijj^in)^ > 0, since the eigenvalues of the collective modes must be real and positive. The 
system becomes instable when the softening is achieved at the critical density rici- 



uj_(T,nci) = 0. 



(25) 
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Taking into account the temperature variation of ujq{T) as explained in Sec. Ill Fl and cal- 
culating self-consistently radius r^. by fl2^ . we find that nci{T) varies from 2 to 6 MPM 
depending on temperature. The result is only different slightly from that obtained in 
(see, Fig. ([5])), since we did not use (HM in [271]. 

Rigourously, such an instability may be interpreted in two ways: either some electrons 
escape from their cavities at nci{T) and that corresponds to the onset of metallization; 
either the system of solvated electrons becomes instable with respect to the formation of a 
state having a finite (non- vanishing) permanent dipole momentum {m.{uj = 0)) 7^ 0. This 



possible phase was first recognized by Turkevich and Cohen 



62 
68 



63] as an excitonic insulator 

(EI), and later studied by different authors [6j, l65|, l66|, l67|, l68|, IfiOj. However, these studies 
treated neutral atoms, which are completely different from the present case (in particular, 
an eventual permanent momentum for the solvated electrons would necessitate a complete 
reorganization of ammonia molecules around the cavities). 

Regarding the onset of metallization within the framework of a generalized GH criterion, 
we calculate the dielectric constant of the system by (l22l) and find that the real part of the 
dielectric constant diverges (zero denominator) when: 



47rnoo(T)ar(0) 

3^00 



(26) 



which defines the second critical density noo{T), i.e. the locus of the polarization catastrophe. 



The similar criterion was considered by Herzfeld [29|], except that the polarizability Oq of 
non-interacting particles was considered. By comparison, using the bare polarizability in 
our formula, we obtain n^Q to be three times larger than that calculated with the use of the 
effective polarizability ar, (we find 14 MPM instead of 5 MPM at T = -70°C). This clearly 
indicates the fundamental role played by the induced dipole-dipole interactions in the MNM 
transition. The comparison between nci{T) and noo{T) shown in Fig.Q indicate that they 
are quite close, i.e. UdiT) ^ UooiT) and corresponds to the onset of metallization. 



C. Dielectric function and optical properties 

As explained above, using Eqn fl22|) we can extract the dielectric function of the solution. 
An example of the calculated e{uj) is presented in Fig. [61 As it is seen, the characteristic 
frequencies range having a nonzero imaginary part and a strong augmentation of the real 
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part of the dielectric function broaden as the concentration rises. At the same time, the 
maximal values of imaginary part ei{uj) and the real part e,.(u;) enormously increase with the 
concentration, and shift to lower frequencies up the the point of polarization catastrophe at 
which er{uj) diverges at = 0. Beyond this point, i.e. for n > rioo ~ 5 MPM, the calculations 
are meaningless since the system is metallic. The calculated real part of the static dielectric 



constant may be also compared with experimental data extracted from the literature [7d. l71|. 
The comparison is shown in Fig. [3 We find a good agreement between the theory and the 
experiments, although slight deviations are observed at high temperatures, probably due to 
the presence of metallic electrons ignored in the present study. 

Using the obtained data on the frequency-dependent dielectric function e{T,uj), we may 
also calculate the optical absorption coefficient A{uj) as: 

A{u;) = Aou;[^ej{u) + eKoo) - e,(a;)]l/^ (27) 

where Aq is a constant of the proportionality, while e.r(co') and ei^uj) are respectively the real 
and imaginary parts of the dielectric function. Figure [8] indicates the calculated concentra- 
tion dependence of A{uj) at the locus of its maximum cjmax at T = — 65°C. This maximum 
decreases at increased metal concentrations, indicating a red shift of the absorption maxi- 
mum observable in the experiments A slight deviation between the theory and the 
experiments is however observed. Thus, this red shift give additional evidence in favor of 
the GH scenario, because it is a signature of the polarization catastrophe |30|. 



D. Thermodynamical properties 

We examine the van der Waals part A/^. At low concentrations the obtained value is close 
to our estimation f|T8l) . although it pronouncedly deviates from it for larger concentrations 
(see Fig. [9]). The first and the second order corrections are significant, which clearly indicates 
the many-body character of the dispersion contribution. We have also calculated all the 
contributions of the total excess free energy, which are depicted in Fig. [101 As it is seen, 
the total change in the free energy is about 10//3 and weakly varies with the concentration: 
the variation does not exceed in the range 0.5MPM < n < 5 MPM. The situation 
is similar with the electrostatic contribution A/c, which only provides an almost constant 
amount of about —2/6^0" to the free energy. At the same time, the van der Waals and 
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the nonpolar contributions (essentially the ideal term in (5)) significatively vary, and their 
delicate balance results in a nonmonotonic behavior of the excess chemical potential ^Utot = 
Afid + AyUc + and the excess pressure Ap, both having their maximum (Fig. [TOl) 
at ns{T) ^ 1.5 MPM for T = — 70°C. The spinodal curve ns{T) of the corresponding 
'liquid-gas' transition is determined by the condition dAfi/dn = 0. Along the curve ns{T), 
the excess compressibility k diverges, indicating the 'liquid-gas' phase separation for larger 
concentration. The influence of the electrostatic interactions to the curve ns{T) is minor, 
they only slightly shift it towards lower concentrations. This is consistent with our basic 
hypothesis: the van der Waals contribution plays the most important role in diluted MAS. 
The spinodal curve ns{T) is represented in Fig. [5], where we have let the van der Waals 
contribution independent on temperature. That gives the locus of the thermodynamical 
instability of the nonmetallic phase. 

At the same time, we may estimate the critical temperature T^^ corresponding to the 
liquid-gas transition with the use of the relations dAfitot(T^^) / dn = d"^ AntotiT^^) / dri^ = 0. 
As a result, we flnd T^^ ~ 500'^K and n^f ~ 11.4 MPM. In fact, the liquid-gas separation 
is hidden by the quantum MIT at rici > ns(T), and the actual state in competition with 
the low density state of localized electrons is not a dense state of solvated electrons, but a 
metallic state of free electrons. Leaving aside details, we have estimated the critical point of 
the phase diagram as the crossing between ns{T) and nci{T). It yields nc{Tc) ~ 2.8 MPM 
and Tc ~ 260 K (see Fig. [5]). Of course, the critical concentration ns{T) signiflcatively 
depends on the size of ions, because the actual packing factor increases with ion size (Fig. 
[TTl) . whereas the critical concentration nc{Tc) varies insignificantly with the changes in the 
ion size. 

Using analytical expressions f|T8l) and (25), we have evaluated the variation of the equilib- 
rium cavity size (see Fig. [T2|) . At low concentrations the obtained value is close to rg = 3.2 
A. The cavity size decreases initially and then pronouncedly rises for metal concentrations 
above 0.5 MPM, the variation in the radius being about 5% at n ^ Q MPM. The excess 
molar volume, which was experimentally measured , is given by Ve = ATc{ai/2 + re)^ in 

our model and the concentration variations of the cavity radius re(n) essentially correspond 
to the the variation of the excess molar volume Ve{n). It is interesting to notice that the ob- 
tained behavior of Ve{n) is experimentally observed for excess molar volumes at temperature 
above Tc . We believe that the initial decrease in may also result from formation 



18 



of eventual multi-electron cavities or due to presence of metallic droplets. 

Finally, we depict the experimental phase diagram of Na-NHs solutions and our critical 
lines nci{T) and ns{T) in Fig. [131 As it is seen these lines give upper estimates with 
respect to the experimental ones, both for the corresponding line of metallization and for 
the phase separation density range. Our calculations indicate that critical concentration 
UgiT) substantially depends on the ion size (Fig. [11]), the phase separation range decreases 
as the size of ion increases due to an enhanced excluded volume effect, while the phase 
separation disappears in the case of Cs^, which is experimentally observed |l|. 



IV. DISCUSSION AND CONCLUSIONS 



Using simple calculations based on a semi-continuum treatment of metal-ammonia solu- 
tions and applying the RPM for interacting excess electrons and ions, we have investigated 
the insulating phase of MAS. Our study has shown that at finite metal concentrations the 
physics of electron solvation is quite different from that of the infinitely diluted case. The 
latter is mainly determined by the cavity formation where the electron is localized, and the 
excess electron can be considered from this point of view as a charge in a cavity surrounded 
by the oriented solvent molecules, while the collective behavior of the electrons is mainly 
controlled by electrostatic interactions. However at intermediate concentrations about sev- 
eral MPM, the excess electrons are governed by dispersion interactions leading to strong 
anomalies in the dielectric response and concentration changes in the absorption maximum. 

Another consequence of dispersion attractions between the electrons, that we have not 
discussed in the paper, is their plausible association. The detailed study of such an asso- 
ciation and related clustering effects is beyond the scope of our current work, because it 
requires an account of the presence of the metallic phase. Nevertheless, our preliminary 
evaluations indicate that a microemulsion phase could arise in the range of concentration 2 
MPM< n < 10 MPM due to the small difference between the free energies of metallic and 
nonmetallic states. This microemulsion phase would be characterized by a large variety of 
aggregates including various stripes, bubbles and so on, as reported in numerical simula- 
tions 19!]. Although magnetic properties of diluted MAS are also beyond the scope of our 
work, we note that the spin-pairing of the localized electrons in the cavities could be a first 
stage of the association indicated above. Does the spin-paring effect result in pure bipolaron 
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formation or in more complicated forms? There is no clear experimental evidence. Some ex- 
perimental data on the spin-pairing were treated as an equilibrium between single electrons 
and spin-pairing bipolarons 73|] , while the concentration changes in magnetic susceptibility 
and excess molar volume were fitted with the use of an assumption of the formation of 
e~ M^e~ complexes |7^]. Quite generally, some deviations of our results with respect to the 
experimental data (on the absorption spectrum and the excess molar volume) may be ex- 
plained in terms of these possible associations. This effect complicates detailed calculations, 
however the formation of such associations does not change our main conclusion about the 
role of induced dipole-dipole interactions, especially for the MNM transition mechanism. 

Our calculations of the dielectric, optical, and thermodynamic properties are confirmed 
by the experimental data, but the most remarkable observation of the dominance of disper- 
sion interactions is the existence of two kinds of instabilities in the nonmetallic phase. The 
first one is a purely quantum transition and related with the onset of metallization, which 
is revealed by a sharp increase in the static dielectric constant when approaching the polar- 
ization catastrophe density. Although we have considered here only the nonmetallic phase, 
we think that the theory of the MNM transition in MAS is actually more complicated and 
provided by a combination of exchange-correlations effects in the metallic phase and dipolar 
interactions between solvated electrons in the insulating phase. This will be the scope of 
our future papers. The second transition is classical and is revealed by a phase separation 
in MAS at low temperatures. This instability is similar to that observed in classical fluids 
with van der Waals attractions. These attractions become dominant at large enough con- 
centrations. At the same time, in our case of MAS the metallic phase is also instable at low 
temperatures for light alkali metals as we have shown in 27|. Therefore, the phase sepa- 
ration in MAS has a dual nature. An additional argument in favor of our dual scenario is 
that both instabilities are to disappear in Cs — NH^ solutions. The latter is experimentally 
confirmed by the absence of the phase separation in these solutions. 

The dominance of dispersion attractions also results in the negativity of dielectric function 
e{uj) over a wide frequency range at /c = 0. At this stage, we have not extended our 



3G, 
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32l |. Together with 



calculations at A; 7^ 0, but e(0, k 7^ 0) should be also negative 
the instability of the metallic phase, that may lead to unusual properties of MAS. For 
example, frustrating the phase separation in fast frozen metal-ammonia solutions, as did 
Ogg 751] a long time ago, the negativity of the dielectric function may result in unusual 
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electronic states with strong attractions [3l| . Is it superconducting or highly conducting but 
nonstationary the answer depends on dynamics and conditions of the electron localization, 
and the existing experimental data on anomalous conductivity in frozen MAS has been 



debated during decades [75|, l76|, l77|. However, this question of possible negativity of the 



dielectric constant has been recently revived by its experimental observation in an expanding 
metal (rubidium) 78| , as well as by the observation of a first-order MNM transition in 
expanded fluid Hg \79l\ . 

Finally, we think that the proposed scenario may be applied to systems involving in- 
teracting quantum particles having a strong polarizability. In our model the behavior of 
quantum polarizable particles is controlled by two dimensionless parameters: ci = ao/rlsoo 
and C2 = Snrf. The first of them is the ratio between the polarizability to the excluded 
volume, while the second one is the relative fraction of the excluded volume in the liquid. 
The parameter ci is small in the case of ordinary polarizable fluids, hence the parameter C2 
should be close to unity to provoke the polarization catastrophe and the present scenario 
remains unlikely for such systems. However in case of self-trapped quantum particles like 
excitons and polarons, Ci may become large enough to provoke a polarization catastrophe 
at low concentrations opening a region in which both the insulating and the metallic phases 
are not stable. The metallization may be associated with a phase separation in that case. 
Hence, our scenario is quite general and may take place in other systems, for instance, in 
alkali metal-alkali halide solutions where solvated electrons, phase separation, and dielectric 



anomalies were experimentally observed [82j. This dual nature could be also the origin of the 
existence of two critical points in the phase diagram of excitons in semiconductors 83|, l84|. 

At this stage the model is not complete, since we have restricted ourselves to low metal 
concentrations. Nevertheless, our calculations have indicated that MAS is an example of 
quantum-classical system whose thermodynamic and dielectric properties are controlled by 
dispersion interactions between self-trapped quantum quasiparticles. We believe that our 
scenario is more general and may be revealed for other quantum particles such as excitons in 
semiconductors and polarons in oxides, which may provide a new insight into metal-insulator 
transitions in condensed matter physics. 
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Appendix 1. Free energy of excess electrons in the insulting phase 

We treat MAS as a system consisting of excess electrons, of monovalent metal ions, 
and A^s classical solvent particles, whose distribution depends not only on their coordinates 
= {Ri, R2, ....Rat^} but also on the orientations w'^^"^ = {wi, W2, ....wjy^} of their 
dipole momenta m. The system of interacting classical particles and electrons is described 
by the grand partition function S given by 

S = {{exp[-f3{H - fi,N - fiiN - ^.AT,)]) J^^. , (28) 

where the symbols (...)gand (..•)g_,_j denote the averages over solvent, electronic, and ionic 
degrees of freedom respectively, H is the total Hamiltonian of the system, fie, fJ-i, and /xare the 
chemical potentials of electrons, ions, and solvent particles. We write the total Hamiltonian 
as the sum of electronic (-f^e), ionic {Hi) and solvent (Hg) contributions: 

NNs 1^1 ^ 

H = T+Y^ Ues{r-^-RsJ, wj) + ^ E | Z 1 + E Mei(r_i-R+,) + (29) 

^ N NNs I Ns 

- E^iidR+i-R+jl) + E Uis(R+i-Rsj,Wj) + -'Y uss(Rsi-'Rsj,Wi - Wj), 

The first term in the right side of (1291) is the kinetic energy of electrons, the second one is 
due to electron-solvent interactions, the next terms are the electrostatic interactions between 
electrons and ions, while the last two terms result from the ion-solvent and solvent-solvent 
contributions, respectively. The electron-solvent term includes electrostatic as well as non- 
electrostatic contributions. The similar contributions are in the ion-ion {uu), the ion-solvent 
(uis), and the solvent-solvent (ugs) potentials. We take into account in (129!) that the inter- 
actions between electrons are screened by the high-frequency dielectric constant £00- 

In principle, the average over quantum and classical degrees of freedom can be provided 
with the use of A^— electron wave function \I'(ri, r2...rAr) or expressed in terms of path inte- 
grals, however neither the computation of the electron wave functions nor the evaluation of 
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path integrals is easy to perform. The most popular version for treating electron-electron 



interactions is based on the local density approximation [85| in which the interactions are 
represented by the sum of electrostatic interactions between classical charges related with 
electron density ne(r)and a short-range contribution due to change in the chemical potential 
/ie(^e(r)) caused by exchange-correlation effects. Unfortunately, such way is not suitable 
for our case, because it ignores the long-range nature of electron-electron correlations. In- 
stead of it, we take into that the electrons are localized and their deviations uj from their 
equilibrium positions R_j = r_j — uj of their centers of mass are small, i.e., uj << R-j. 
Then expanding the electron-electron and electron-ion interactions with respect to this small 
parameter and restrict ourselves by the dipole approximation, we obtain 

The first terms in these relations are the Coulomb interactions between classical charges 
(cations and anions), the second terms are responsible for interactions between permanent 
dipoles and classical charges, while the last term in (130|) is due to interactions between in- 
duced dipoles, but the last term in fl3Tl) is interactions between classical charges and induced 
dipoles. The contribution responsible for interactions between the permanent dipoles and 
classical charges vanishes after averaging over quantum degrees of freedom, because the lo- 
calized electrons have spherical symmetry in the ground state and their average momenta 
(m) oc (u) are equal to zero. The similar situation for the contribution caused by the 
last term in fl3T]) . it disappears after averaging over the classical degrees of freedom due to 
isotropic distribution of classical charges. As a result, expansions (130!) and (13T!) ) provide 
only the classical electrostatic interactions (the first terms in (l30l)-( l3T|) ) and the induced 
dipolar interactions (the last term in (l30l) ). 

The set {Ri^\ R^^\ R|^=^, wf^"-^} corresponds to classical degrees of freedom. Av- 
eraging over the classical degrees, we obtain an effective hamiltonian H^f depending on 
correlation functions between solvent and excess charges and on quantum coordinates. The 
evaluation of the correlation functions can be performed by the theory of integral equations 
of classical liquids [86| . Such approach has been apphed in [48|, |49|] to model a mixture of 



dipolar solvent particles, excess electrons, and ions. However, these equations can be solved 
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only numerically, and to simplify the calculations, we apply the RPM and treat the solvation 
of excess charges as a formation of cavities and screening of charges localized in the cavities. 
In this case the effective hamiltonian Hej is written as 

1 ^ 

Hef = Nho + AF,i + — ^ u, ■ T(R_, - R_,) • u,-, (32) 

where ho is the one-electron contribution for noninteracting electrons, AF^i is the change 
in the free energy due to the classical interactions between solvated charges (cations and 
anions), while the last term corresponds to dispersion interactions between the electrons. 
The average of the one-electron hamiltonian Hq over quantum degrees of freedom can be 



obtained by the integral equation methods for an excess electron in a polar liquid 20|, l87| , 



and the free energy /o can be expressed [2l|] as: 



/o = -^/^o(r)V2M/o(r)rfr+^^ / \'^li^)Kf{r - r,)\^liT,)\drdr, + U,, (33) 

where \&o(r) is the electron wave function in the ground state, fcav is the nonelectrostatic 
contribution due to cavity formation, while Me/(r) is the effective interaction potential due to 



polarization of solvent molecules. It can be found 471] that the long-range part of the effective 



potential has the asymptotic Uef{r — > oo) = —(e^— e ^) /r, while cavity contribution can be 



parameterized as a sum of volume and surface terms [21[ . Finally, using the semi-continuum 
approximation for the effective potential Me/(r), we obtain ([2]), while the change AF/N in 
the total free energy is written by ([1]). 

Appendix 2. The effective polarizability of a quantum polarizable fluid 

Let us consider a large sample of M >> 1 Drude oscillators in a volume V with surface S 
and extract a sphere containing particles from this large sample (with M >> N » 1). 
In the spirit of the Clausius-Mossotti calculation, we may evaluate the polarization of the 
independent sphere subjected to an external field, which is nothing but the Lorentz local 
field fl2T|) . The momentum induced by the external field on i-th particle is: 

TV 

= J2[ao'l-e^'T,,]-' ■ Ee{uj), (34) 
j 

while the expression for the average moment {m{uj)) is written as 

(m(^)) = /f^K 'l-£^'T,,r'\ ■ M^) = a{u)Ee{u)/e^, (35) 
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where (...) is the average over distributions of particles, while the last identity in the above 
relation is the definition of renormalized effective polarizability a{uj) in which we account 
that the external field is screened by the solvent containing in the sample. Then the change 
AP in the macroscopic polarization, determined as AP =Y,f t^ii^^), is written in terms of 
effective polarizability by ( l20l) . that allows us to obtain (1221) . 

The averaging procedure results in the self-consistent equation for the renormalized po- 

3e evaluated by the integral equations method 
3, Q provides the analytical expression (fT2|) 



arizability written formally as f[TT|) . It can 



88l . l89| , while the Pade approximation 
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for this function. Substituting it into ( ITTl) . we obtain quadratic equation ( |T3|) for a{u) whose 
solution is complex, while their real and imaginary parts are given by: 

ar{uj) = aoB{uj)A-\u)[^Jl - 2A{uj)B-^{uj) - 1], (36) 

ai{uj^ <u<uj+)= aoA'\uj)^j2A{uj) - B^{uo), (37) 

where B{uj) = lu'^luq'^ — l + baos^ , A{ij) = 2aQ£^b[uj'^tUQ'^ — l + ab^^ao£^], and the functions 

a and b are determined in terms of dimensionless parameter x = Snr^ as: 

TTx J - O.3I682; - 0.3205x2 + 0.1078a;^ Trn ^ , , 
= 2¥|f (1 - 0.5236^)2 ] = 3^1'°^^^' (^^^ 

(2.70797+ 1.68918X -0.3157x2) 57m/i(x) 

b(x) = 0.72505n^ , ! = ^V- (39) 

^ ' (1 - 0.59056X + 0.20059x3)/o(x) 8/o(x) ^ ^ 

The root of (1371) is equal zero at edge eigen- frequencies c<j_and u;+ : 

ulin) = ul + b{n)e-^ ± 2a^'\n)e-^, (40) 

determining the low and the high edge frequencies of the DOS spectrum. Hence, we can 
calculate effective polarizability a{(jj) and the DOS D{uj) for a given concentration n and 
frequency cuq, while the later can be derived from rg by ([3]). 

The solution of f|T3|) depends on the two dimensionless parameters Ci = I/EooT^^j^q and 
C2 = 8nrg. To obtain analytical estimations of the expressions, we consider the limiting 
case C2 — i> 0, or more precisely rg — > (called as the point oscillator limit). We have 
a = Sirn J^^r~^dr = Tm/Sr^ and b = 57in/8, which leads to expressions (fTSl) and f|T6|) for 
the imaginary part of the effective polarizability and the edge frequencies. At the same time, 
the real part ar{0) becomes in the same limit: 

a,(0,n^O)=ao[l+ „ ^ 3 ]' (41) 
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while it is equal to 0:^(0, rici) = 2ao/(l ~ Svrnciao/S^oo) at the critical concentration rid. 
The instability occurring at rid is given by uj^{nci) = 0, which leads to: 

~ 4n ~ (2re)3 ' ^^^^ 
where Ci(T) = {6/7i){eooC^{T)/re)'^ is the numerical factor weakly depending on tempera- 
ture. This factor Ci is about of 1/3 at low temperatures, and hence the critical concentration 
is roughly rid ~ l/3(2re)'^. Substituting the analytical expressions for rid and ar{0,nd) we 
may evaluate the second instability occurring at riao given by: 

noo _ ao(l - 57rncooao/8£oo) ^ 
rid 2eoor^ 

The two instabilities in this limit are found to be quite close to that of metallization, since 
noo/ud ~ 0.92 at low temperatures. 

Now we evaluate the zero-point energy of the collective mode in the limiting case. Sub- 
stituting (fTSll and (flGjl into the expression of the zero-point vibrations, we have 

A/. = ^ [" - u'f/^u^ - ulf'du - ^. (44) 

Introducing the new variable which is related with the frequency u) as 

cu^ = cUq + cos ?/, (45) 

where 5 = (47m/3£:^rg)^/^, we can rewrite the integral as 

A/, = ^ r[(l + 6u,'cosyy/' - 1] sin' ydy. (46) 

TT Jo 

Expanding the root in the integrand into the series with respect to parameter ScOq'^, i.e. 
(1 + ScOq'^ cosyY^'^ — 1 = Sujq'^ cosy / 2 — {5uJq'^ cos?/)^/8 + we obtain 



Expressing 5 via density n, we have the final expression ( 1181) for the dispersion contribution. 
Completing the consideration of the limiting case, we indicate that the above estimate for 
the energy of zero-point vibrations is only asymptotically valid at C2 — ^ 0, however we should 
account the higher contributions at finite n, moreover the series is divergent at the point rid 
of the polarization catastrophe, since SoJq'^ = 1. 
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FIG. 1: Schematic sketch of modcL The solvated cations witli charge +e are represented by the 
orange circles and have cartesian coordinates R+j of their centers of mass, whereas the cavities 
with localized electrons are represented by the blue circles and have the coordinates R_i of their 
centers of mass. The instantaneous position of i—th electron is r_i = R_i + Ui, where Ui is the 
relative coordinate of the electron with respect to the cavity center. Both cations and solvated 
electrons are assumed to have the same size a. The inset represents the solution in which the 
solvent is treated as a continuum medium with dielectric constant eNHzi^)- 
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FIG. 2: The dimensionless DOS D{io,n)ujQ versus dimensionless frequency uj/lvq at various metal 
concentration n, T = — 70*^0, and fixed ojq = 0.9 eV and re = 3.2A. The values of metal concen- 
tration n are indicated at the corresponding lines. 




FIG. 3: The dimensionless imaginary part ai{uj) / ao{0) (a) and the real part ar{uj)/ao{0) (b) of 
effective polarizabihty at various metal concentration n, T = — 70'^C, and fixed ujq = 0.9 eV and 
Te = 3.2A. The notions are the same as in Fig. [2j 
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FIG. 4: The dimensionless edge square frequencies uPl/ujq (solid lines) and uj\/ojq (dashed lines) 
versus metal concentration at fixed cavity radius tq (thick lines) and with the account of concen- 
tration dependence of the radius re(n) (thin lines). The notions are the same as in Fig. [5J 




FIG. 5: Critical lines for Na-NHa solutions. The solid lines correspond to the instability concentra- 
tion nci(T), the dashed ones to the spinodal line n.s(T) (see text), the dotted line to divergent static 
dielectric constant (polarization catastrophe). The thick solid line indicates nci(T) calculated at 
fixed cavity radius = 3.2 A and frequency ujq{T), while the thin one to that with the account of 
the concentration dependencies of the above parameters. 
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FIG. 6: The real part €r{uj) (a) and the imaginary part ei{uj) (b) of the dielectric function at various 
metal concentration n, T = — TO'^C, and fixed = 0.9 eV and = 3.2A. The notions are the 
same as in Fig. [2j 
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FIG. 7: The static dielectric constant er{oo = 0) versus the metal concentration; the circles symbols 
correspond to the experimental data [toI on the low-frequency dielectric constant at T = -|-20''C, 
the squares to that at T = — 35°C [71]. The dashed and solid curves show our results at T = -|-20°C 
and —35°C, respectively. The thin and the thick curves correspond to calculations at fixed cavity 
radius rg and with the account of the concentration dependence of the radius re(n), respectively. 
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FIG. 8: The locus of the maximum t^max of optical absorption in MAS, the triangle symbols indicate 
the experimental data on the absorption maximum in Na-NHa at T = — 65°C obtained from [t^] 
and square symbols from [73], while the solid curve shows our results at the same temperature. 
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FIG. 9: Concentration dependence of the dispersion contribution Affi{n) calculated at T = —70*^0. 
The dashed line corresponds to the contribution normalized with respect to its limiting values with 
account the first-order correction A/^i = —TrnuJoaQ/We'^r^, while the solid line to that with 
respect to the second-order correction A/rf2 = + 5vrnao/24e^rg). All other conditions are 

the same as that in Fig. [2l 
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FIG. 10: The dimensionless contributions to the excess free energy (a), to the excess chemical 
potential (b) at T = —70*^0. The nonpolar contributions are depicted by solid lines, the dispersion 
ones by the dashed lines, the electrostatic ones by the dotted lines, and the total change by the 
dashed-dotted lines. The arrow indicates the critical concentration at which the compressibility 
is equal to infinity. All other conditions are the same as that in Fig. [2j 





FIG. 11: Influence of the ion radius on the locus of the spino dal p oint n^. The symbols corresponds 
to the van der Waals radii of the relevant ions, derived from 80|. All other conditions corresponds 
to that in Fig. [H 




FIG. 12: Concentration dependence of the cavity radius rg{n) at various temperatures. 



42 



10 



-10 



-30 



-50 



-70 



T(C) 



4 




• 


exp [23] 




cubic 




"cl(T) 




"3(1) 


T 


exp [24] 




exp[1] 



I / 
I • 



n(MPM) 



8 



FIG. 13: Phase diagram of Na-NHs solution. The ex per imental data on the locus of the phase 
separation are indicated by squares 



24l ]. respectively, while the dashed-dotted 



231 ] and triangles 

line indicates the cubic interpolation of these data. The diamonds show the change in sign of the 
derivative of the conductivity coefficient da/dT, which is used to estimate the locus of the MNM 
transition The dashed curve corresponds to the low bond of the spinodal. The solid curve 
depicts our calculations of the locus of the polarization catastrophe. The bottom of the figure 
corresponds to the solidification of ammonia. 
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